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To exploit Rashba effect in a 2D electron gas on silicon surface for spin transport, it is necessary to have 
surface reconstruction with spin-split metallic surface-state bands. However, metals with strong spin-orbit 
coupling (e.g., Bi, Tl, Sb, Pt) induce reconstructions on silicon with almost exclusively spin-split insulating 
bands. We propose a strategy to create spin-split metallic bands using a dense 2D alloy layer containing a 
metal with strong spin-orbit coupling and another metal to modifj^ the surface reconstruction. Here we 
report two examples, i.e., alloying Bi/Si(lll)\/3 x \/3 reconstruction with Na and Tl/Si(lll)l X 1 
reconstruction with Pb. The strategy provides a new paradigm for creating metallic surface state bands with 
various spin textures on silicon and therefore enhances the possibility to integrate fascinating and promising 
capabilities of spintronics with current semiconductor technology. 



The Rashba spin splitting^ in the two-dimensional electron gas systems on semiconductors is considered to 
be the key concept for many promising spintronics applications. To combine Rashba-effect based spin- 
tronics with a silicon technology metallic surface states with a strong spin-orbit coupling on a silicon surface^'^ 
are demanded. Unfortunately, the most of the metal/silicon systems with spin- split surface state bands (e.g., 
Bi/Si(lll)A/3xA/3^-^Tl/Si(lll)l X 1 Sb/Si(lll)A/3 x and Pt/Si(110)"6" X 5 dive semiconducting 
The only known exceptions are the Tl/Si(lll)l X 1 modified by adsorption of additional TP^ and Au/Si(lll) 
X modified by adsorption of In, Tl, Cs or Na^. It is worth noting that in both cases adsorbates do not alter 
the basic atomic arrangement of the pristine surface and their effect resides just in eliminating or generating 
surface defects (as for Au/Si(lll) and Tl/Si(lll), respectively) and doping electrons to the available surface-state 
bands. In the present study, we propose a novel more universal strategy for tailoring the spin-split metallic surface 
states. The main concept resides in formation of the 2D alloy layers containing a metal with a strong spin-orbit 
coupling and another suitable metal to obtain dense reconstructions with spin-split metallic bands. Validity of the 
approach is demonstrated by the formation of the Bi-Na and Tl-Pb 2D alloys on Si(lll) surface possessing the 
required properties. 

Results 

Figure 1 (left panel) summarizes the main structural and electronic properties of the pristine Bi/Si(lll) 
X surface. It contains 1.0 ML (monolayer, 7.8 X 10^^ cm"^) of Bi arranged into the trimers centered at the 
T4 sites on the bulk- truncated Si(lll) surface^^. The main features of its electronic structure are spin-split non- 
metallic bands, SI (SI') and S2 (S2')^. Room-temperature (RT) deposition of 0.33 ML Na onto this surface results 
in the formation of the ordered Bi-Na alloy layer whose structural and electronic properties are shown in Figure 1 
(right panel). Comparison of its atomic structure with that of the original Bi/Si(lll)A/3xV^ surface demon- 
strates that the Bi trimers are preserved but they increase in size by —8% (Bi-Bi bond length changes from 3.13 to 
3.34 A) and twisted by ±6.8°. The resultant arrangement of Bi atoms at the surface can be thought then as a 
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Figure 1 | Structural and electronic properties of the Bi/Si(ll l)\/3 x ^3 (left panel) and (Bi, Na) /Si(lll)\/3 x \/3 (right panel), (a) and (e) 30 X 30 

STM images. Upper and lower halves in (e) present experimental and simulated empty-state ( + 0.3 V) STM images, respectively, (b) and (f) Structural 
models where Bi atoms are shown by violet circles, Na atoms by green circles, and Si atoms by gray circles. Positions of the outlined \/3 x \/3 unit cells in 
(e) and (f) correspond to each other. The inset in (b) shows a sketch of reciprocal space geometry with boundaries of the \/3 x \/3 and 1X1 SBZs given by 
solid black and dashed blue lines, respectively, (c) and (g) Experimental ARPES spectra and (d) and (h) calculated electronic band structures. In the 
calculated band structures, the size of green circles corresponds to the strength of the surface character summed over all orbitals at a particular h value. 



honeycomb network of the chained Bi trimers with Na atoms occu- 
pying the T4 sites in the center of each honeycomb unit. 

Changes in atomic structure of the surface cause modification of 
its electronic band structure as revealed by angle-resolved photoelec- 
tron spectroscopy (ARPES) measurements and density-functional 
theory (DFT) calculations (Figure Ig and h). We would like to 
remark that a perfect agreement between the calculated band struc- 
ture and that determined in the experiment as well as coincidence of 
the experimental and simulated STM images can serve an indication 
of the proper (Bi, Na)/Si(lll) structural model. The most essential 
feature in the Bi-Na alloy band structure is the appeared metallic 
spin-split surface-state band denoted by SI (ST). 

Consider now the Pb-induced modification oftheTl/Si(l 11)1 X 1 
surface. Figure 2 (left panel) presents the main structural and elec- 
tronic properties of the pristine Tl/Si(lll)l X 1 reconstruction. It 
contains 1.0 ML of Tl atoms occupying every T4 site on the bulk- 
truncated Si(lll) surface^^"^^. The band structure of the Tl/Si(lll)l 
X 1 surface includes the insulating SI surface- state band whose 
occurrence was proved before both experimentally^'^^ and theoret- 
ically^'^^. Spin -splitting of this band has been concluded from the 
spin-resolved ARPES results^ and from previous^ and current theor- 
etical calculations (Figure 2d) and is clearly resolved in the spin- 
unpolarized ARPES spectrum in Figure 2c. Additional feature can 
be noted, namely, the shallow metalUc band S2 around the iC point. It 
has recently been recognized that electron filling of this band is 
associated with extra Tl atoms on the Tl/Si(lll)l X 1 surface^^ in 
the form of specific surface defects 



Adding 0.33 ML of Pb to the Tl/Si(lll)l X 1 surface at RT pro- 
duces a homogeneous Tl-Pb alloy having x \/3 periodicity. Its 
atomic arrangement (Figure 2f) looks qualitatively akin that of the 
Bi-Na alloy (Figure If), namely the first adsorbate, Bi or Tl, form the 
chained trimers arranged into a honeycomb network, while atoms of 
the second adsorbate, Na or Pb, occupy the centers of the honeycomb 
units. The difference is that Pb atom residence and Tl trimers centers 
in the Tl-Pb alloy are in the Ti (on-top) sites while the positions of Na 
atoms and Bi trimers in the Bi-Na alloy are the T4 sites. Both Bi-Na 
and Tl-Pb alloys are confined essentially in single atomic layers 
with Na being by —0.62 A higher than Bi and Pb by 0.37 A higher 
than Tl. 

Substantial structural changes due to Tl/Pb alloying, including, in 
particular, the change in the surface periodicity, 1 X 1 to \/3 x \/3, 
leads to changing the band structure. The most essential newly 
developed features are the two spin-split metallic surface-state bands 
(denoted by YA (ST) and S2 (S2')) which are clearly seen in 
Figure 2g and h. 

The calculated Fermi surfaces of the alloys reproduce nicely the 
contours observed experimentally (Figure 3). The main features of 
the alloy spin texture predicted by DFT calculations are also illu- 
strated in Figure 3. The planar spin components show in -plane heh- 
city, suggestive of a Rashba scenario. For momentum vectors along 
the Y—K directions, the spin is fully aligned in-plane and perpen- 
dicular to the momentum vector. The out-of-plane spin component 
undulates between positive and negative values along the contours 
according to the Cs^ symmetry of the surfaces. 
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Figure 2 | Structural and electronic properties of the Tl/Si(lll)l X 1 (left panel) and (Tl, Pb) ^ Si(lll)\/3 x (right panel), (a) and (e) 30 X 30 A" 

empty-state ( + 0.4 V) STM images. Upper and lower halves in (e) present experimental and simulated empty-state ( + 0.4 V) STM images, respectively, 
(b) and (f) Structural models where Tl atoms are shown by violet circles, Pb atoms by green circles, and Si atoms by gray circles. Positions of the outlined 
V3 X ^/3 unit cells in (e) and (f) correspond to each other, (c) and (g) Experimental ARPES spectra and (d) and (h) calculated electronic band structures. 



Si(111)V3xV3-Bi,Na 

(a)i.o- 





l(b> 








A / t 


O +Sz 

o -Sz 




Si(111)V3xV3-TI,Pb 

(e)i.o- 








\ t 
\ ^ 


^ / V 

O +Sz 
O -Sz 








Figure 3 | Calculated band structure of the (Bi,Na)/Si(lll)A/3 x a/3 (left panel) and (Tl,Pb) y^Si(lll)\/3 x \/3 (right panel) including spin-orbit 
coupling illustrating the (a, e) in-plane and (b, f) out-of-plane polarization components, (c, g) Experimental and (d, h) calculated Fermi contours. Inset in 
(g) shows the fine structure within the outlined area around the M point with a greater contrast. Arrows adjacent to the calculated contours and their 
length denote the in-plane spin component. The out-of-plane spin component is indicated by the colour with red and blue corresponding to the upward 
and downward directions, respectively. White colour indicates fully in-plane spin alignment. The tic values in parentheses in (d) refer to the pocket 
contours around the K. 
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Figure4| Calculated electron density maps for (a) Bi/Si( 111) \/3 x \/3, (b) (Bi, Na)/Si(lll)\/3 x v^, (c) Tl/Si( 111)1 X 1, and (d) (Tl, Pb)/Si(lll) 
a/3 X a/S surfaces. The multi-colour (blue-yellow) scale for electron density is used to highlight the interconnected electron density network at metallic 
alloy surfaces. 



In particular, the (Bi, Na)/Si(lll) band structure contains spin- 
split El (ET) band which shows up in the Fermi map (Figure 3c and 
d) as two hexagon- shaped contours with corners pointing in the 
T—M directions (where the splitting of the band is maximal, 
momentum splitting AA;|| = 0.044 A"^ and energy splitting AE^ = 
210 meV). Besides, the SI (ST) band forms hole pockets having the 
shape of smoothed triangles around the K points. The out-of-plane 
spin component for pocket contours is negligible (albeit finite and 
different for neighboring pockets), the in-plane component domi- 
nates and the counterclockwise spin helicity is the same for all six 
pocket contours. 

As for the (Tl, Pb)/Si(lll), it displays the two spin-split metallic 
bands, SI (SI') and i:2 (2:2'). In the Fermi map of the SI (SI'), the 
outer contour has almost round shape, while the inner contour is a 
hexagon with corners pointing in the T—K directions. The maximal 
splitting for the SI (SI'), Ak^ = 0.038 A"' and AEp = 250 meV, is 
along T—M direction. The other spin-split metallic band S2 (S2') 
shows up as hexagonal contours with corners along the T—M dir- 
ection, i.e. rotated by 30° compared to the inner contour of the SI 
band. The maximal splitting for the S2 (S2'), AA;|| = 0.050 A"' and 
AEp =140 meV, is along T—K direction. 

Comparing the Fermi contours of the present alloys with those 
found for spin-split surface states of the reported phases with hexa- 
gonal symmetry, one can notice certain common features. For exam- 
ple, the contour with two hexagons having corners in the T—M 
directions (as for the (Bi, Na)/Si(lll) SI (SI') and (Tl, Pb)/Si(lll) 
S2 (S2') bands) was also observed at the Pb/Ge(lll)A/3 x \/3 sur- 
faced The contour of (Tl, Pb)/Si(lll) SI (SI') band with hexagon 
having corners in the T—K directions inside a circle is akin that re- 
ported for Au/Si(lll)y3 X 73^ and Au/Ge(lll)y3 x ^3 

Discussion 

The shown examples of 2D alloying provide a hint for elucidating the 
pathways for converting the originally semiconducting surface into 
the metallic one. Note that both original surfaces, Bi/Si(lll) 
a/3 X a/3 and Tl/Si(lll)l X 1, contain 1.0 ML of metal adsorbate. 
Taking into account that the surface lattice constant of Si(l 1 1) is 3.84 
A, such a density of adsorbate atoms is not sufficient to ensure the 
overlapping of their electron wave functions. This is clearly seen in 
the calculated electron density maps in Figure 4a and c: electron 
density is localized either around Bi trimers or single Tl atoms, 
respectively. Note that it is typical for metal-induced Si(lll) recon- 
structions where the atomic layers with metallic properties develop 
usually at metal coverage beyond 1.0 ML (e.g., In/Si(lll)A/7 x a/3 
with 1 .2 and 2.4 ML W or Pb/Si( 1 1 1 ) a/3 x a/3 with 1.33 ML Pb"'). 
Another possibility is realized when the metal layer of 1.0 ML is 
compressed into the interconnected network due to Si atoms incor- 
porated within the same reconstructed layer (e.g., Ag/Si(lll)A/3 x 
a/3 or Au/Si(lll)A/3 X a/3 "^'"^ each containing 1.0 ML metal 



and 1.0 ML Si). One can see that adding 1/3 ML Na to the 
Bi/Si(lll)A/3 X a/3 produces a similar effect: Na adsorption causes 
originally isolated Bi trimers to rotate and grow in size until they 
form a chained-trimer structure. As a result, electron density assoc- 
iated with Bi atoms becomes arranged into the interconnected net- 
work (Figure 4b) responsible for appearance of the surface metallic 
properties. The (Tl, Pb)/Si(lll) system looks seemingly similar, as 
originally isolated Tl atoms also form chained-trimer structure. 
However, the electron density map (Figure 4d) shows that actually 
the Tl- associated electron density is interconnected through that of 
the Pb atoms. Hence, appearance of the surface metallic properties is 
controlled here just by increasing metal-atom coverage at the surface. 

Thus, our results demonstrate that formation of 2D alloys contain- 
ing a metal with a strong spin-orbit coupling allows to obtain spin- 
split metallic surface-state bands on silicon. However, finding an 
appropriate metal pair is actually a challenging task. For example, 
while Na is a suitable component for making 2D alloy with Bi, it 
appears that the other alkali metals, Li and Cs, are not alloyed with Bi 
and grow as 3D islands on ^5 — Bi/Si(lll)A/3 x a/3 surface. Never- 
theless, the list of the prospective alloys displaying variety of spin- 
related properties is believed to be very vast. 

Methods 

Sample preparation. The STM and ARPES experiments were performed in an ultra- 
high- vacuum Omicron MULTIPROBE system with a base pressure better than —2.0 
X 10"^° Torr. Atomically- clean Si(lll)7 X 7 surfaces were prepared in situ by 
flashing to 1280°C after the samples were first outgassed at 600°C for several hours. 
Pristine Bi/Si(lll)\/3 x a/3 and Tl/Si(lll)l X 1 reconstructions were formed by 
depositing 1 ML of the corresponding species onto Si(lll)7 X 7 surface held at 
~450°C and ~300°C, respectively. Sodium was deposited from a well-outgassed 
commercial SAES chromate dispenser, lead from a heated Pb-stuffed Mo tube. 

STM. STM images were acquired using Omicron variable-temperature STM-XA 
operating in a constant- current mode. Electrochemically- etched W tips and 
mechanically cut Ptir tips were used as STM probes after annealing in vacuum. 

ARPES. ARPES measurements were conducted in the ultrahigh vacuum chamber 
Omicron MULTIPROBE using VG Scienta R3000 electron analyzer and high-flux He 
discharge lamp {hv = 21.2 eV) with toroidal-grating monochromator as a light 
source. 

Ab initio random structure search. Atomic structure of the new alloyed 
(Bi, Na) y/si( 1 1 1 ) + and (Tl, Pb) /si( 1 1 1 ) \/3 x \/3 reconstructions were 
elucidated using ab initio random structure search technique^*^ after it had been 
successfully tested for the known Bi/Si( 111) a/3 X a/3 andTl/Si(lll)l X 1 structures. 

First-principles calculations. Our calculations were based on DFT as implemented 
in the Vienna ab initio simulation package VASP^^'^*^, using a planewave basis and the 
projector- augmented wave approach^'' for describing the electron-ion interaction. 
The generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof 
(PBE)^° has been used for the exchange correlation (XC) potential. The Hamiltonian 
contains the scalar relativistic corrections, and the spin-orbit interaction (SOI) was 
taken into account by the second variation method as has been implemented in VASP 
by Kresse and Lebacq^\ To simulate the reconstructions we use a slab consisting of 12 
bilayers (BL). Hydrogen atoms were used to passivate the Si dangling bonds at the 
bottom of the slab. Both bulk Si lattice constant and the atomic positions within the 
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three BLs of the slab were optimized including SOI self- consistently. The silicon 
atoms of deeper layers were kept fixed at the bulk crystalline positions. 
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